AFRL- VS-HA-TR-2005-1 162 


Observations  and  Modeling  of  the  Upper  Mesosphere: 
Mesopause  Characteristics,  Inversion  Layers,  and  Bores 


Peter  P.  Wintersteiner 
Edward  Cohen 


ARCON  Corporation 
260  Bear  Hill  Road 
Waltham,  MA  02451-1080 


Scientific  Report  No.  1 
26  October  2005 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


AIR  FORCE  RESEARCH  LABORATORY 
Space  Vehicles  Directorate 
29  Randolph  Road 

Am  FORCE  MATERIEL  COMMAND 
Hanscom  AFB,  MA  01731-3010 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 
AFRL-VS-HA-TR-2005-1 162 


/signed/  /signed/ 

JEREMY  WINICK  ROBERT  BELAND,  Chief 

Contract  Manager  Battlespace  Surveillance  Innovation  Center 


This  report  has  been  reviewed  by  the  ESC  Public  Affairs  Office  (PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS). 

Qualified  requestors  may  obtain  additional  copies  from  the  Defense  Technical 
Information  Center  (DTIC).  All  others  should  apply  to  the  National  Technical 
Information  Service. 

If  your  address  has  changed,  if  you  wish  to  be  removed  from  the  mailing  list,  or  if 
the  addressee  is  no  longer  employed  by  your  organization,  please  notify 
AFRL/VSIM,  29  Randolph  Rd.,  Hanscom  AFB,  MA  01731-3010.  This  will  assist 
us  in  maintaining  a  current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a 
specific  document  require  that  it  be  returned. 

Using  Government  drawings,  specifications,  or  other  data  included  in  this  document  for 
any  purpose  other  than  Government  procurement  does  not  in  any  way  obligate  the  U  S. 
Government.  The  fact  that  the  Government  formulated  or  supplied  the  drawings, 
specifications,  or  other  data  does  not  license  the  holder  or  any  other  person  or 
corporation;  or  convey  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented 
invention  that  may  relate  to  them. 

This  report  is  published  in  the  interest  of  scientific  and  technical  information  exchange 
and  its  publication  does  not  constitute  the  Government’s  approval  or  disapproval  of  its 
ideas  or  findings. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  (or  this  collection  o(  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  inducing  suggestions  for  redudnQ 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway.  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  It  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY) 
10-26-2005 

2.  REPORT  TYPE 

Scientific  Report  No.  1 

3.  DATES  COVERED  (From  -  To) 

08-01-2004  to  07-31-2005 

4.  TITLE  AND  SUBTITLE 

Observations  and  Modeline  of  the  Upper  Mesosphere:  Mesopause  Characteristics.  Inversion 

5a.  CONTRACT  NUMBER 

FA87 1 8-04-C-003 1 

Layers,  and  Bores 

5b.  GRANT  NUMBER 

N/A 

5c.  PROGRAM  ELEMENT  NUMBER 

612301 

6.  AUTHOR(S) 

Peter  P.  Wintersteiner  and  Edward  Cohen 

5d.  PROJECT  NUMBER 

2301 

5e.  TASK  NUMBER 

BD 

51.  WORK  UNIT  NUMBER 

AI 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

AJR.CON  Corporation 

260  Bear  Hill  Road 

Waltham,  MA  02451-1080 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESSES) 

Air  Force  Research  Laboratory 

29  Randolph  Rd. 

Hanscom  AFB,  MA  01731-3010 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AFRL/VSBYM 

11.  SPONSOR/MONITOR  S  REPORT 
NUMBER(S) 

AFRL-VS-H  A-TR-2005- 1 162 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Using  data  from  TIMED/SABER,  we  studied  the  global  behavior  of  the  mesopause  at  solstice  (northern-hemisphere 
summer).  The  existence  of  the  two-level  structure  with  an  abrupt  altitude  transition  near  25N  is  confirmed.  Its  altitude  in  the 
“summer”  zone  is  shown  to  increase  with  increasing  (northerly)  latitude;  in  the  “winter”  zone,  it  is  found  somewhat  below  the 
accepted  100-km  level  at  midlatitudes.  Temperatures  are  in  the  165-1 80K  range  throughout  the  winter  zone  but  north  of  25N 
they  fall  continuously,  reaching  -120K  at  80N.  We  show  that  mesospheric  temperature  inversion  layers  (TILs)  are  very 
common,  and  that  large  TILs  affect  the  mesopause  characteristics.  We  also  characterize  the  occurrence  distribution  of  TILs  in 
latitude  and  local  time. 

We  studied  the  properties  of  mesospheric  bores  by  numerically  solving  the  differential  equations  governing  their  propa¬ 
gation  modes.  Using  an  impulsive  forcing  function  and  measured  wind  and  temperature  profiles,  we  determined  the  vertical 
displacement  as  a  function  of  time  and  horizontal  distance  for  different  altitudes  within  the  duct,  and  thereby  were  able  to 
study  the  evolution  of  the  wave.  We  were  able  to  reproduce  the  phase  speed  of  an  observed  bore  to  a  good  approximation. 

15.  SUBJECT  TERMS 

Atmospheric  radiance  structure;  mesospheric  temperature  inversion  layer;  MIL;  TIL;  SABER;  global  mesopause;  gravity  waves’  ducting; 
mesospheric  bores;  Navier  Stokes  equations 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Dr.  Jeremy  Winick 

a.  REPORT 

UNCLASSIFIED 

b.  ABSTRACT 

UNCLASSIFIED 

c.  THIS  PAGE 

UNCLASSIFIED 

UU 

53 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

781-377-3619 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


Table  of  Contents 


1 .  Introduction  1 

2.  Investigations  of  the  MLT  Region  Using  SABER  Data  1 

2. 1 .  Global  Characteristics  of  the  Mesopause  at  Solstice  2 

2.2.  Mesospheric  Temperature  Inversion  Layers  6 

3.  Modeling  Mesospheric  Bores  9 

3.1.  Introduction  9 

3.2.  Formulation  of  Taylor-Goldstein  and  Helmholtz  Equations  10 

3.3.  Altitude-Limited  Data  Considerations  1 1 

3.4.  Solution  of  Taylor-Goldstein  and  Helmholtz  Equations  12 

3.5.  The  Benjamin-Davis-Ono-Burgers  (BDOB)  Equation  and  its  Solution  14 

3.6.  Vertical  Parcel  Displacements  1 7 

3.7.  Summary  of  Results  17 

References  19 


iii 


List  of  Figures 

1.  SABER  mesopause  altitude  distribution  as  a  function  of  latitude  for  21  June-9  July  2003.  22 

2.  Similar  to  Figure  1,  but  using  all  ascending  data  from  5  July  2003  only.  22 

3.  Distribution  of  nighttime  mesopause  altitudes,  by  latitude,  in  late  June  and  early  July  2003.  23 

4.  SABER  mesopause  altitude  distribution  as  a  function  of  latitude  for  5-17  July  2003.  24 

5.  Temperature  profiles  for  three  successive  scans  on  3  July  2003.  24 

6.  SABER  mesopause  temperature  distribution  for  the  conditions  of  Figure  1 .  25 

7.  SABER  mesopause  temperature  distribution  for  the  conditions  of  Figure  2.  25 

8a.  Mean  temperature  profiles  for  ascending  data  on  5  July  2003,  5-degree  latitude  bins  35-75°N.  26 

8b.  As  in  Figure  8a,  but  for  latitude  bins  35°N-5°S.  27 

8c.  As  in  Figure  8a,  but  for  latitude  bins  5-45°S.  28 

9.  Comparison  of  mean  temperature  within  five  degrees  of  the  equator  on  27  March  2002.  29 

10.  Four  nighttime  temperature  profiles  from  27  March  2002.  30 

1 1 .  Four  nighttime  temperature  profiles  from  13  July  2002,  from  consecutive  events.  30 

12.  TIL  occurrence  probabilities,  by  latitude  and  local  time,  for  two  yaw  cycles  in  2002.  3 1 

13.  Mean  altitude  distribution  ofTILs,  by  latitude  and  local  time,  for  the  yaw  cycles  in  Figure  12.  31 

14.  TIL  occurrence  probability,  as  in  Figure  12  but  for  one  yaw  cycle  at  solstice.  32 

15.  Mean  TIL  altitudes,  plotted  against  local  time,  for  equatorial  events  during  March-May  2002.  32 

16.  Helmholtz  equation  propagation  constant  (squared).  33 

17.  Helmholtz  equation  eigenmode  #2  and  phase  speed,  from  the  rational  Chebyshev  polynomial  34 

solution. 

18.  Convergence  of  phase  speed  estimates  using  rational  Chebyshev  polynomials.  35 

19.  Solution  A(X,T)  of  the  BDOB  equation.  36 

20.  Snapshot  of  the  solution  A( X,  T  =  0,100)  of  the  BDOB  equation.  37 

21.  Time-staggered  snapshots  of  the  solution  A(X,T)  of  the  BDOB  equation.  38 

22.  Snapshot  of  the  solution  A (X,T  =  0,100)  of  the  BDOB  equation.  39 

23.  Solution  A(X,T)  of  the  BDOB  equation  assuming  a  very  long  initial  pulse.  40 

24.  Snapshot  of  the  solution  A(X,  T  =  0, 100)  of  the  BDOB  equation.  41 

25.  Evolving  vertical  displacements  at  various  altitudes,  from  the  Taylor-Goldstein  and  BDOB  42 

equations. 

26.  Evolving  vertical  displacements  at  various  altitudes,  from  the  Taylor-Goldstein  and  BDOB  43 

equations. 

27.  Evolving  vertical  displacements  at  various  altitudes,  from  the  Taylor-Goldstein  and  BDOB  44 

equations. 

28.  Evolving  vertical  displacements  at  various  altitudes,  from  the  Taylor-Goldstein  and  BDOB  45 

equations. 


v 


1.  INTRODUCTION 

This  Interim  Scientific  Report  describes  work  that  was  undertaken  during  the  first 
twelve  months  that  contract  #FA8718-04-C-0031  was  in  effect.  The  period  covered  is 
from  1  August  2004  through  31  July  2005.  The  topics  we  discuss  fall  into  two  categories, 
both  pertaining  to  the  characteristics  and  dynamics  of  the  upper  mesosphere  and  their 
effects  on  the  optical  emissions  that  arise  in  that  region. 

The  first  category  involves  data  from  the  Sounding  of  the  Atmosphere  using  Broad¬ 
band  Emission  Radiometry  (SABER)  instrument  [Russell  et  al„  1999 ]  on  the  TIMED 
satellite.  SABER  is  a  ten-channel  radiometer  that  is  sensitive  to  emissions  from  six 
important  molecular  constituents  in  the  infrared  (IR)  and  near-infrared  (NIR)  wave¬ 
lengths.  The  measured  radiance  profiles  are  used  to  retrieve  temperature  and  pressure, 
volume  emission  rates  (VERs),  and  constituent  mixing  ratios.  These  and  other  products 
have  been  used  in  many  studies  involving  the  energy  balance  of  the  mesosphere  and 
lower  thermosphere  (MLT),  the  target  region  for  the  TIMED  mission.  One  subsection  of 
this  report  deals  with  global  characteristics  of  the  mesopause  region  at  solstice,  as 
revealed  by  these  data.  Mesospheric  temperature  inversion  layers  (TILs)  are  treated  in  the 
other  subsection.  Mesospheric  TILs  have  been  observed  at  various  locations  in  the  past, 
but  SABER  provides  the  first  global  perspective  on  this  subject,  and  reveals  that  they  are 
often  very  large  and  much  more  common  than  was  previously  realized. 

The  second  category  involves  numerical  solutions  of  equations  that  govern  the  propa¬ 
gation  of  gravity  waves  within  ducts  in  the  upper  mesosphere.  In  some  circumstances, 
TILs  provide  such  ducts,  and  thereby  support  the  propagation  of  large  amplitude  solitary 
waves  or  internal  bores,  structures  that  are  analogous  to  bores  that  have  been  observed  in 
the  troposphere  and  in  tidal  estuaries  (rivers).  Such  structures  propagate  over  great  hori¬ 
zontal  distances  and,  in  turn,  produce  wavelike  optical  structure  that  has  been  seen  in 
several  emissions  of  the  mesopause  region.  In  addition  to  basic  understanding  of  how 
they  propagate,  the  goals  are  to  quantify  the  likely  range  of  the  temperature  and  radiance 
excursions  that  might  be  seen  and  the  rapidity  with  which  observed  radiance  can  be 
expected  to  vary. 

2.  INVESTIGATIONS  OF  THE  MLT  REGION  USING  SABER  DATA 

SABER  [Russell  et  al.,  1999]  is  a  ten-channel  radiometer,  one  of  four  instruments 
aboard  the  TIMED  satellite.  Since  early  in  2002,  it  has  been  continuously  scanning  the 
limb  between  the  ground  and  tangent  heights  of  approximately  300  km,  recording  IR  and 
NIR  emissions  from  CO2,  H2O,  O3,  NO,  OH,  and  O2. 

The  SABER  CO2  15  pm  limb  radiance  is  used  to  determine  atmospheric  temperature 
and  pressure  up  to  ~  1 00  km.  This  is  much  higher  than  was  previously  possible,  and  is 
made  possible  by  a  new  non-LTE  retrieval  algorithm  [Mertens  et  al.,  2001;  2002;  2004]. 
Because  of  good  S/N  ratios  in  the  15  pm  channels,  the  algorithm  is  actually  stable  up  to 
-125  km,  but  the  specified  accuracy  of  three  Kelvins  in  temperature  is  not  achievable  at 
those  altitudes  because  several  quantities  are  uncertain.  Minor-constituent  densities,  cool¬ 
ing  rates,  and  volume  emission  rates  (VERs)  are  also  returned  for  various  altitude  ranges 
in  the  middle  atmosphere.  The  radiance  data  are  referred  to  as  Level  IB  data;  tempera¬ 
ture/pressure,  VERs  and  mixing  ratios  are  Level  2  data. 
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Because  of  the  continuous  duty  cycle,  the  high  inclination  of  the  TIMED  orbit,  and  the 
sensitivity  of  the  detectors,  there  is  a  prodigious  flow  of  data,  and  near-global  coverage  is 
assured  for  many  products  up  to  quite  high  altitudes.  Our  studies,  reported  below,  take 
advantage  of  this  broad  coverage  to  examine  seasonal  and  latitudinal  variations  of  impor¬ 
tant  characteristics  of  the  mesosphere  and  lower  thermosphere  (MLT).  Among  these  are 
the  principal  mesopause  parameters — its  altitude  and  temperature.  The  global  behavior  of 
these  parameters  at  solstice  is  the  subject  of  the  first  subsection  (2.1.).  We  are  also 
studying  mesospheric  temperature  inversion  layers  (TILs),  which  are  now  recognized  as 
relatively  common  occurrences  everywhere  except  in  the  high-latitude  summer  region. 
This  is  reported  in  the  second  subsection  (2.2.). 

In  both  of  these  cases,  as  in  many  other  aspects  of  work  on  the  MLT  region,  atmos¬ 
pheric  tides  comprise  an  important  influence  that  must  be  accounted  for.  Therefore,  a 
considerable  effort  has  been  made  to  distinguish  the  data  we  use  according  to  local  time 
(LT).  We  often  speak  of  a  “parameter  space”  consisting  of  season,  latitude,  and  local 
time,  and  we  have  devoted  some  effort  to  extracting  and  sorting  the  retrieved  temperature 
and  measured  radiance  accordingly.  In  this  context,  it  is  important  to  note  that  the  pre¬ 
cession  of  the  TIMED  orbit  is  such  that  local  times  of  the  equator-crossings  decrease  by 
about  .2  hours  per  day.  It  takes  almost  60  days  to  view  all  accessible  local  times,  between 
the  ascending  and  descending  portions  of  the  orbit.  The  ascending/descending  LT  differ¬ 
ence  at  the  equator  is  about  8.8  hours.  Also,  times  near  local  noon  can  never  be  observed 
because  of  detector  sensitivity.  Therefore,  whenever  the  satellite  precesses  so  the  instru¬ 
ments  point  to  within  a  certain  angle  close  to  the  sun,  a  yaw  maneuver  is  performed  and  a 
new  ~60-day  “yaw  cycle”  is  begun. 

2.1.  Global  Characteristics  of  the  Mesopause  at  Solstice 

We  report  here  on  a  preliminary  study  of  the  global  variability  of  mesopause  height 
and  temperature  using  SABER  Level  2  (version  1 .04)  temperatures.  This  effort  was  con¬ 
current  with  our  investigation  of  MLT  temperature-inversion  layers,  described  below 
[also,  Picard  et  al.,  2004;  Winick  et  al.,  2004],  The  layer-locating  algorithm  we  devel¬ 
oped  for  that  purpose  required  finding  the  mesopause  as  a  first  step.  We  have  concentra¬ 
ted  on  the  northward-looking  yaw  cycle  of  May-July  2003  that  includes  the  solstice 
(northern  hemisphere  summer)  period.  Eventually,  this  will  be  extended  to  a  full  year’s 
worth  of  data. 

For  the  current  study,  the  mesopause  has  been  defined  as  the  location  of  the  coldest 
temperature  above  the  stratopause.  Figure  1  shows  the  distribution  of  nighttime  meso¬ 
pause  altitudes  as  a  function  of  latitude  for  all  SABER  events  between  21  June  and  9  July 
with  local  time  between  1  and  2  hours.  The  results  show  the  low  “summer”  mesopause  at 
northern  high  and  mid-latitudes,  the  high  “regular”  mesopause  extending  through  the 
winter  hemisphere,  and  an  abrupt  transition  between  the  two  at  ~25°N,  in  accord  with 
previous  observations.  Figure  2  shows  the  altitude  distribution  using  all  of  the  “night¬ 
time”  data  from  5  July  only.  Superposed  upon  the  individual  event  data  in  Figure  2  are 
the  mesopause  altitudes  calculated  from  mean  profiles  using  all  nighttime  events — those 
near  local  midnight — in  5-degree  latitude  bins.  Taking  altitudes  from  the  mean  profiles 
emphasizes  the  sharpness  of  the  transition.  The  mean  altitude  falls  between  94  and  101 
km  everywhere  south  of  the  transition. 
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These  results  are  generally  in  accord  with  past  observations,  for  example  those  made 
during  the  Polarstern  cruise  of  1996  [von  Zahn  et  al.,  1996 ].  Those  nightly-averaged 
measurements  led  to  the  notion  of  the  global  bipolar  mesopause,  at  ~86  km  in  the  mid-  to 
high-latitude  summer  and  ~100  km  everywhere  else,  that  is  now  well  accepted  [e.g., 
Thulasiraman  and  Nee,  2002].  A  followup  GCM  model  study  [ Berger  and  von  Zahn, 
1999]  demonstrated  that  the  two-level  mesopause  is  an  inherent  feature  of  the  atmos¬ 
phere,  not  an  artifact  of  nighttime-only  measurements,  and  elucidated  the  roles  of  certain 
physical  mechanisms  (IR  cooling,  chemical  heating,  momentum  deposition  by  gravity 
waves)  that  are  crucial  for  producing  it.  It  also  showed  that  tides  affect  the  mesopause 
altitude  and  temperature,  but  do  not  themselves  cause  the  two-level  phenomenon. 

Although  the  principal  characteristics  of  the  SABER  altitude  distributions  agree  with 
these  and  other  previous  measurements,  a  more  detailed  picture  emerges  as  a  result  of  the 
extensive  coverage  and  the  large  number  of  separate  events  that  are  available  for  study. 
For  one  thing,  the  spread  in  altitudes  can  be  determined  globally.  Based  on  results  in  our 
figures,  for  instance,  this  spread  is  greater  than  the  ~3  km  originally  suggested  by  von 
Zahn  et  al.  [1996],  especially  south  of  the  transition  point.  This  is  true  even  if  one 
disregards  those  events  for  which  the  mesopause  appears  in  the  -75-85  km  range;  those 
may  represent  anomalous  circumstances,  but  they  are  common  enough  to  require  an 
explanation.  The  “regular”  mesopause  also  appears  to  be  lower  by  almost  5  km  at  sub¬ 
tropical  latitudes  in  the  winter  hemisphere  than  at  the  equator  and  at  midlatitudes.  Finally, 
the  mean  altitudes  shown  in  Figure  2,  on  both  sides  of  the  transition,  are  lower  by  at  least 
2  km  than  the  commonly  accepted  altitudes  of  the  “summer”  and  “regular”  mesopause. 

The  spread  of  ±3  km  noted  above  [vo/i  Zahn  et  al.,  1996]  represents  a  variation  of  dif¬ 
ferent  nights’  averages  at  changing  latitudes.  It  can  be  compared  to  the  variation  of  the 
SABER  mean  altitudes  given  by  the  red  squares  in  Figure  2,  which  is  also  approximately 
±3  km.  Those  points,  however,  are  centered  about  97  km  rather  than  100  km.  They  also 
appear  to  have  a  systematic  dependence  on  latitude  rather  than  an  approximately  constant 
value,  producing  a  lower  (~94  km)  mean  mesopause  near  ~5-25°S.  This  pattern  is 
repeated  in  other  subsets  of  the  May-July  data  for  local  times  around  midnight,  but  not 
for  times  much  earlier  in  the  evening  or  later  in  the  morning.  Note,  these  mean  altitudes 
were  each  determined  from  a  single  average  temperature  profile  for  events  within  a  parti¬ 
cular  LT/latitude  bin,  not  by  averaging  the  individual  events’  mesopause  altitudes.  There¬ 
fore,  the  fact  that  the  distributions  shown  in  Figures  1  and  2  are  also  skewed  to  lower  alti¬ 
tudes  near  ~5-25°S  is  a  separate  indication  that  some  physical  mechanism  lowers  the 
mesopause  there;  it  is  not  a  redundant  confirmation  of  the  pattern  seen  in  the  mean 
altitude. 

Figure  3  gives  a  different,  and  perhaps  more  useful,  view  of  the  spread  of  mesopause 
altitudes.  For  6-degree  bins  spaced  throughout  the  range  of  latitudes,  it  shows  the  distri¬ 
bution  for  events  near  local  midnight  during  the  latter  part  of  June  and  early  July.  The  last 
six  frames  depict  the  high  “regular”  mesopause  exclusively,  and  one  can  see  a  spread  of 
~5  km  about  the  centroid  in  each  case.  In  contrast,  the  first  four  frames  show  the  low 
“summer”  mesopause.  At  the  most  northerly  latitudes,  the  distribution  is  very  narrow,  no 
more  than  ~2  km  about  the  centroid,  but  it  broadens  as  the  latitude  approaches  the  transi¬ 
tion  region. 
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It  is  particularly  noteworthy  that,  in  the  summer  high  latitudes,  the  mesopause  altitude 
shown  in  Figures  1  and  2  increases  by  about  5  km  between  40^  and  82°N.  This  was 
hinted  at  in  the  Polarstern  data  and  predicted  by  the  aforementioned  model  study  [Berger 
and  von  Zahn,  1999].  However,  this  is  the  first  time  that  such  a  result  has  been  confirmed 
at  far  northern  latitudes,  past  69°. 

The  polar  summer  mesopause  altitudes  shown  here,  and  appearing  in  other  subsets  of 
SABER  data,  are  consistently  lower  by  ~2-3  km  than  is  typically  observed  [ Lubken  et  al., 
1999].  This  is  easy  to  see  because  of  the  narrow  distribution  there,  but  from  Figure  1 — 
and  more  clearly  from  other  SABER  data — the  regular  mesopause  also  appears  to  be  ~2- 
3  km  below  its  nominal  height  of  100  km.  As  a  result  of  this,  the  question  of  whether 
there  is  a  bias  in  the  SABER  altitudes  has  been  raised.  A  review  of  the  pressure  registra¬ 
tion  procedure,  which  could  have  induced  an  error  of  that  magnitude,  was  conducted,  but 
the  changes  that  were  instituted  appeared  to  have  a  minimal  effect.  To  verify  this,  we 
compared  two  days’  data  processed  with  the  preliminary  version  1.06  software  with  some 
of  the  older  vl.04  data.  The  new  results  reflect  a  number  of  changes  in  addition  to  the 
pressure  registration,  but  the  new  mesopause  altitudes  were  within  ~2  km  of  the  old  ones 
for  almost  all  events.  More  importantly,  there  was  no  systematic  difference;  the  distribu¬ 
tions  were  practically  unchanged.  Therefore,  the  discrepancy  between  SABER  meso¬ 
pause  altitudes  and  previous  measurements,  particularly  in  the  summer  polar  region, 
remains  a  mystery. 

In  Figure  1,  a  small  but  not  negligible  number  of  events  in  the  winter  hemisphere  have 
very  low  mesopause  altitudes.  In  similar  format,  Figure  4  shows  the  altitude  distribution 
in  the  afternoon  (LT  15-16  h)  using  data  taken  during  July.  In  this  case,  a  much  larger 
number  of  events  south  of  the  transition  zone  have  this  characteristic  and,  as  a  result,  the 
transition  zone  itself  is  far  less  distinct.  In  contrast,  during  morning  hours  (7-9  LT;  not 
shown)  there  are  few  very  low-mesopause  events.  This  apparent  dependence  on  local 
time  suggests  a  tidal  influence.  The  tides  are  particularly  strong  in  the  upper  mesosphere 
and  they  perturb  the  temperature  profiles  considerably.  Part  of  our  future  effort  will  be  to 
model  these  effects  explicitly. 

In  order  to  investigate  this,  we  looked  at  many  individual  events  throughout  the  yaw 
cycle.  We  found  a  lot  of  events  with  what  appears  to  be  strong  tidal  or  other  wave  struc¬ 
ture  superposed  upon  what  one  might  suppose  to  be  the  mean  state.  (This  is  especially  so 
in  the  winter  hemisphere.)  This  raises  the  question  of  whether  the  definition  of  the  meso¬ 
pause  as  the  super-stratopause  temperature  minimum  is  reasonable.  Three  retrieved  temp¬ 
erature  profiles  shown  in  Figure  5  illustrate  the  problem.  They  come  from  three  consecu¬ 
tive  events  in  the  southern  hemisphere,  from  orbit  8492  on  July  5.  They  each  have  two 
temperature  inversions  and  three  minima  between  the  stratopause  and  the  thermosphere. 
Despite  their  apparent  similarity,  we  derive  very  different  mesopause  altitudes  of  90,  67, 
and  103  km  for  them.  For  event  83  in  particular,  the  one  with  the  highest  mesopause,  the 
three  widely-separated  minima  are  less  than  one  degree  different  in  temperature,  as  indi¬ 
cated  in  Figure  5.  Many  other  events  have  structure  quite  like  these,  and  some  of  them 
produce  the  low-mesopause  points  that  we  find  in  the  altitude  plots. 

The  question  still  remains  as  to  what  produces  large  temperature  inversions  like  those 
in  Figure  5.  Tidal  amplitudes  are  not  as  large  as  the  TILs  shown  here,  but  nonlinear 
wave-wave  interactions  involving  the  tides  are  thought  to  be  capable  of  augmenting 
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them.  Whatever  the  exact  process  is,  it  seems  clear  from  these  and  many  other  SABER 
profiles  that  wave  effects  have  a  profound  effect  on  the  retrieved  temperature  and  hence 
on  the  mesopause  parameters  derived  from  them.  The  existence  of  a  very  low-altitude 
mesopause  is  most  likely  due  to  such  wave  activity.  This  argues  for  a  redefinition  that 
requires  averaging  not  only  over  short  periods  to  eliminate  gravity  wave  perturbations, 
but  also  over  longer  periods  as  well  to  eliminate  the  effect  of  tides  [ Berger  and  von  Zahn, 
1999].  Given  the  slow  precession  of  the  TIMED  orbit,  of  course,  the  latter  procedure 
would  drastically  compromise  the  SABER  temporal  resolution. 

Figures  6  and  7  show  the  mesopause  temperatures  corresponding  to  the  two  sets  of 
events  depicted  in  Figures  1  and  2.  The  cold  polar  summer  mesopause  is  apparent  in 
each,  as  is  a  gradual  (rather  than  abrupt)  transition  from  “summer”  to  “regular”  condi¬ 
tions.  The  full  range  of  temperatures  in  the  winter  hemisphere  is  quite  large,  -40  K.  The 
temperature  of  the  mean  profiles  varies  by  -15  K  over  the  latitude  range  20°N-50°S  for 
the  July  5  data.  The  spread  of  temperatures  north  of  the  transition  region  is  about  half  of 
what  it  is  in  the  south. 

Figures  8a-c,  using  nighttime  data  from  5  July,  illustrate  how  the  temperature  structure 
changes  through  the  transition  region  between  the  “summer”  and  “regular”  regions. 
These  figures  show  that  the  cold  summer  polar  mesopause  (Figure  8a)  evolves  into  a 
double-minimum  structure  at  northern  mid-  and  low  latitudes  (Figure  8b)  before  transi¬ 
tioning  to  the  high  regular  mesopause  (Figure  8b-c).  Wave  structure  is  noticeable  in  some 
of  the  plots,  despite  the  fact  that  each  average  is  derived  from  20  or  more  individual 
profiles  from  the  full  range  of  longitudes.  This  indicates  that  there  is  similar  structure 
from  event  to  event,  and  it  is  most  likely  driven  by  processes  that  are  consistent  from 
place  to  place— rather  than  sporadic  in  nature  as,  for  example,  gravity  waves  are. 

Studying  characteristics  of  the  mesopause  is  an  obvious  and  important  application  of 
the  SABER  database.  The  major  advance  that  it  represents  is  in  its  near-global  coverage, 
and  the  sheer  number  of  individual  measurements  it  provides  (-1400  events  per  day).  For 
example,  evidence  for  the  existence  of  the  global  two-level  mesopause  has  been  pieced 
together  from  many  past  measurements  made  at  various  locales  by  different  techniques 
and  at  different  times.  SABER  now  provides  a  single  database  capable  of  validating  the 
notion.  The  results  that  we  have  examined  from  this  one  60-day  period  certainly  support 
the  idea.  They  also  show  that,  besides  describing  the  mean  state  of  the  atmosphere  at  a 
given  time  and  locale,  we  have  an  unprecedented  view  of  the  variability  of  the  MET 
region  as  well.  One  finds  that  the  altitude  of  the  mesopause  is  quite  variable,  more  so  in 
the  winter  hemisphere  than  in  the  summer,  but  the  mean  altitude  is  fairly  stable.  One 
finds  an  abrupt  altitude  transition  from  “summer”  to  “regular”  conditions,  but  a  smooth 
transition  for  the  temperature.  In  fact,  SABER  data  provide  an  ideal  tool  for  studying  the 
transition  region,  as  suggested  by  Figure  8.  There  is  a  large  range  of  mesopause  temper¬ 
atures  at  any  given  latitude  in  the  winter  hemisphere,  but  a  much  smaller  range  in  the 
polar  summer.  Moreover,  there  are  a  lot  of  events  displaying  evidence  of  severe  wave 
activity,  which  can  cause  the  mesopause  to  appear  at  anomalously  low  altitudes. 

We  have  also  verified  the  prediction  that  the  altitude  of  the  low  polar  summer  meso¬ 
pause  increases  smoothly  by  -5  km  between  40  and  80°N. 
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The  SABER  data  prove  that  a  complete  picture  of  mesopause  characteristics  requires 
more  than  a  snapshot  in  time  or  space.  Despite  the  rich  and  detailed  information  it  pro¬ 
vides,  however,  the  difficulty  of  distinguishing  between  diurnal  and  seasonal  changes 
remains.  This  problem  arises  because  of  the  slow  precession  rate,  which  allows  us  to 
view  most  local  times  only  once  every  yaw  cycle,  and  the  tides,  which  strongly  modulate 
conditions  of  the  MLT  region  on  an  hourly  basis.  We  think  that  the  best  way  to  resolve 
these  difficulties  is  to  produce  a  detailed  (and  even  redundant,  using  different  years)  view 
of  data  throughout  the  full  parameter  space  (latitude/'LT/season),  and  complement  that 
with  an  effort  to  model  tidal  effects. 

2.2.  Mesospheric  Temperature  Inversion  Layers 

As  noted  above,  the  SABER  database  provides  an  unprecedented  opportunity  to  study 
the  global  behavior  of  mesospheric  temperature  structure.  One  particularly  interesting 
feature  is  the  temperature  inversion  layer  (TIL).  Mesospheric  TILs,  which  have  been  the 
subject  of  great  research  interest  in  recent  years  [e.g.,  Leblanc  and  Hauchecorne,  1997; 
Meriwether  and  Gardner,  2000;  Meriwether  and  Gerrard,  2004;  Sassi  et  al.,  2002; 
Huang  et  al.,  1998;  2002],  appear  with  regularity  in  the  SABER  data.  Each  of  the  temp¬ 
erature  profiles  shown  in  Figure  5  in  the  preceding  section  has  more  than  one  TIL.  As 
noted  there,  wave  structure  may  perturb  the  mesosphere  in  many  ways,  even  to  such  an 
extent  that  the  mesopause  appears  at  very  low  altitudes.  We  have  come  upon  a  number  of 
events  with  TILs  that  are  more  than  twice  as  large  as  the  ones  shown  in  Figure  5,  having 
amplitudes  as  great  as  1 00  K — the  “amplitude”  being  defined  as  the  temperature  differ¬ 
ence  from  an  underlying  minimum  to  a  local  maximum  -5-10  km  higher  up.  This  section 
describes  an  investigation  of  the  TIL  phenomenon  as  seen  in  the  SABER  data. 

The  first  step  we  have  taken  is  to  determine  when  and  where  TILs  are  most  likely  to 
appear.  That  is,  we  wished  to  determine,  for  different  seasons  of  the  year,  the  latitudes, 
local  times,  and  altitudes  where  they  are  seen  and,  conversely,  where  they  are  not  seen. 
The  discrimination  according  to  LT  is  important  because  of  the  role  of  the  tides,  which  is 
not  completely  understood  [e.g.,  Meriwether  and  Gerrard,  2004 ]  but  clearly  a  factor  in 
producing  them. 

This  point  is  illustrated  in  Figure  9.  In  that  figure,  we  juxtapose  the  mean  temperature 
profiles  at  the  equator  for  all  events  on  the  ascending  and  descending  portions  of  the 
orbits  for  a  single  day  in  March  2002,  and  plot  their  difference  as  well.  The  local  times 
for  these  events  were  approximately  0.1  h  and  9  h,  respectively.  The  profiles  are  offset  by 
more  than  50  K  at  two  places,  much  more  than  the  rms  deviations  that  are  also  plotted. 
The  mesopause  altitudes  and  temperatures  are  distinctly  different,  and  a  TIL  of  more  than 
60  K  amplitude  appears  on  the  ascending  (nighttime)  side.  The  difference  plot  reveals  the 
tidal  influence  most  strikingly.  Oscillations  with  a  vertical  wavelength  of  -22-25  km  and 
an  amplitude  increasing  with  altitude  point  to  the  diurnal  tide  as  the  predominant  influ¬ 
ence — as  expected  at  the  equator  at  this  time  of  year — despite  the  fact  that  the  local  times 
are  less  than  12  h  apart.  The  size  of  the  temperature  differences  is  surprising,  since  the 
diurnal  amplitude  is  expected  to  be  only  -10-15  K  in  the  upper  mesosphere  in  March 
[GSWM,  2004]. 

Figures  10  and  1 1  depict  several  nighttime  events  with  TILs.  In  Figure  10,  events  of 
March  27,  2002,  include  one  extremely  large  TIL  (event  4).  One  can  see  that  the  altitudes 
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at  which  these  layers  occur  vary  considerably.  Many  of  the  largest  ones  appear  between 
~83  and  ~93  km  (as  shown  in  Figure  10),  but  the  full  altitude  range  extends  at  least  10 
km  above  and  below  those  points.  We  find  that,  when  TILs  appear  in  event  after  event,  as 
on  the  two  orbits  chosen  here,  the  variation  of  the  altitudes  is  usually  systematic  rather 
than  random.  That  is,  the  heights  change  gradually  with  latitude,  and  frequently  repeat  the 
pattern  on  consecutive  orbits.  Figure  11  shows  four  consecutive  events  with  TILs  all 
centered  near  ~80-84  km,  from  an  orbit  of  July  13,  2002. 

In  order  to  examine  TIL  occurrence  rates,  we  analyzed  data  from  three  complete  yaw 
cycles  in  2002  and  2003.  Two  cycles  began  near  equinox,  one  looking  north  and  one 
south,  and  one  near  solstice  was  north-looking.  There  were  significant  data  drop-outs  in 
each  period,  but  we  were  still  able  to  cover  most  of  the  latitude/LT  parameter  space, 
which  we  separated  into  3-degree  latitude  bins  and  1-hr  LT  bins.  We  determined  the 
percentage  of  events  with  TILs  larger  than  certain  thresholds,  and  the  mean  altitude  of 
such  layers,  for  each  bin  and  for  each  yaw  cycle.  The  thresholds  were  1 0,  20,  and  40  K. 

Probability  results  for  the  two  equinox  cycles  are  given  in  Figure  12  for  20-K  TILs.  It 
shows  that  inversion  layers  are  strongly  correlated  with  LT,  and  also  that  the  variation 
with  LT  depends  on  latitude.  There  appear  to  be  more  TILs  during  the  March-May  period 
than  during  the  September-November  period,  but  in  both  cases  we  find  a  heavy  concen¬ 
tration  throughout  most  of  the  nighttime  near  the  equator.  We  also  find  a  somewhat  light¬ 
er  concentration  in  the  early  morning  and  a  much  lighter  concentration  in  the  early  after¬ 
noon  near  30  degrees  in  the  spring  hemisphere,  plus  some  in  the  early  morning  near  30 
degrees  in  the  fall  hemisphere.  There  are  very  few  at  latitudes  polewards  of  60  degrees. 

One  remarkable  discovery  is  that  TILs  of  this  size  appear  with  a  likelihood  exceeding 
90%  near  the  equator  at  certain  times.  Even  for  40-K  TILs  (not  shown),  it  approaches 
50%  in  the  late  evening. 

Figure  13  gives  the  mean  TIL  altitude  for  these  two  cycles,  in  the  same  format.  For 
most  of  the  bins,  the  distribution  of  altitudes  has  a  single  peak  with  a  half-width  (rms 
deviation)  of  only  ~2  km.  However,  there  are  some  regions  of  the  parameter  space  where 
the  distributions  are  either  quite  broad  or  bipolar,  which  makes  the  quantity  less  meaning¬ 
ful.  Nevertheless,  one  sees  that  the  altitudes  definitely  correlate  with  LT  and  latitude  in 
much  the  same  way  as  the  probabilities  do.  At  the  equator,  for  example,  the  TILs  tend  to 
be  just  above  80  km  at  night  (as  shown  in  Figure  10)  and  much  higher  in  the  daytime 
(where  they  are  also  less  numerous).  At  mid-latitudes,  the  average  nighttime  TIL  is  much 
higher,  also  in  accord  with  Figure  10. 

These  two  yaw  cycles  are  consistent  in  that  there  is  a  fairly  clear  indication  of  diurnal 
patterns  at  the  equator  where  the  diurnal  tide  is  strongest  (especially  near  equinox),  and  a 
somewhat  less  clear  indication  of  semidiurnal  behavior  near  mid-latitudes  where  the 
semidiurnal  tide  tends  to  dominate.  This  can  be  seen  in  both  the  probabilities  and  the 
altitudes  of  the  TILs.  Note  that  the  patterns  north  and  south  of  the  equator  following  the 
March  equinox  are  similar  to  those  south  and  north  of  it,  respectively,  following  the  Sep¬ 
tember  equinox. 

For  the  solstice  cycle  that  we  examined,  May-July  2003,  there  were  numerous  missing 
data.  However,  one  can  see  from  Figure  14  that  the  likelihood  of  finding  TILs  as  large  as 
20  K  is  substantially  less  than  near  equinox  for  all  regions  of  the  parameter  space.  There 
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is  still  a  small  concentration  at  night  at  the  equator,  and  also  in  early  morning  and  late 
afternoon  at  midlatitudes  in  the  winter  hemisphere.  The  mean  altitudes  (not  shown)  of  the 
TILs  that  are  seen  do  vary  with  LT,  but  the  data  are  too  sparse  to  make  much  of  the 
patterns.  There  are  virtually  no  large  TILs  in  the  summer  polar  region,  which  seems  con¬ 
sistent  with  the  fact  that  the  lapse  rate  is  so  large  there. 

Besides  quantifying  the  distributions,  a  second  task  was  to  determine  whether  or  not 
coherent  inversion  structures  commonly  persist  over  extended  distances  (hundreds  or 
thousands  of  kilometers)  in  the  horizontal  direction.  This  is  important  for  several  reasons, 
among  which  are  the  fact  that  a  large-scale  TIL  may  comprise  a  stable  duct  for  long¬ 
distance  propagation  of  gravity  waves,  and/or  support  mesospheric  undular  bores  [ Dewan 
and  Picard,  1998;  2001 ].  This  may  also  shed  some  light  on  the  role  of  gravity  waves  in 
producing  TILs  [ Huang  et  al.,  1998;  Liu  et  al.,  2000].  We  therefore  sought  layers  that 
were  similar  in  amplitude  and  altitude  within  sequences  of  consecutive  events.  Our 
criterion  was  that  qualifying  events  should  all  have  TILs  exceeding  the  chosen  threshold 
(usually  20  K)  at  altitudes  differing  by  no  more  than  3  km. 

The  along-track  direction  of  motion  is  mainly  north-south  at  latitudes  where  most 
TILs  occur,  and  events  are  separated  by  300-400  km.  We  find  that  sequences  of  six  or 
more  20-K  TILs,  representing  -2000  km  or  more  in  distance,  are  very  common,  and  we 
found  several  extending  beyond  4000  km.  Even  with  the  more  restrictive  40-K  require¬ 
ment,  sequences  of  -2000-3000  km  are  common  near  the  equator.  In  general,  the  preval¬ 
ence  of  extended  TILs  mirrors  the  probability  that  they  occur  at  all,  insofar  as  they  are 
much  less  common  or  absent  at  high  latitudes  and  solstice. 

From  this,  one  concludes  that  inversion  layers  very  often  do  appear  as  extensive  struc¬ 
tures,  at  least  in  the  north-south  direction  for  some  latitudes  and  times  of  day.  One  might 
assume  that  the  same  is  true  in  all  directions,  but  it  is  difficult  to  demonstrate  this  because 
consecutive  orbits  are  separated  by  more  than  2000  km  (at  low  latitudes)  and  90  minutes. 
Also,  LT  remains  nearly  unchanged.  Therefore,  even  if  similar  structure  were  observed 
within  a  latitude  region  during  consecutive  events  of  two  or  more  successive  orbits,  one 
could  not  tell  whether  it  was  static  and  long-lived  or  changing  with  the  phase  of  the  tide. 

Understanding  the  role  of  the  atmospheric  tides  in  producing  TILs  is  a  third,  and  very 
important,  reason  for  studying  their  global  behavior.  We  have  begun  a  modeling  effort 
using  the  state-of-the-art  Global  Scale  Wave  Model  [GSWM,  2004 ]  to  try  to  explain  the 
SABER  occurrence  probabilities  and  amplitudes  discussed  earlier.  Although  this  effort  is 
far  from  complete,  we  have  made  some  simple  comparisons  of  TIL  altitudes  with  GSWM 
predictions.  In  Figure  15,  we  plot  centroids  of  TIL  altitude  distributions  against  local  time 
for  events  near  the  equator  during  the  March-May  yaw  cycle.  We  also  plot  the  altitude  of 
maxima  in  the  tidal  temperature  perturbation  from  GSWM.  The  model  altitudes  descend 
regularly  by  about  25  km  in  24  h,  as  expected  for  a  dominant  diurnal  component.  From 
-4  h  to  24  h  LT,  the  TIL  altitudes  also  descend,  somewhat  more  gradually  but  in  approx¬ 
imate  agreement  with  the  model.  However,  from  -0  h  to  8  h,  TIL  altitudes  between  -80 
and  85  km  do  not  vary  with  LT.  (From  4  to  8  h,  there  is  a  bipolar  distribution.) 

This  behavior  appears  to  be  consistent  with  conclusions  drawn  by  Meriwether  and 
Gerrard  [2004].  They  distinguish  two  types  of  mesospheric  TILs:  a  “higher”  TIL,  in  the 
altitude  range  -85-100  km  that  is  tidally  driven  and  amplified  by  nonlinear  gravity-wave 
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interactions;  and  a  “lower”  type,  below  ~80  km,  that  is  caused  by  dissipating  planetary 
waves.  The  “higher”  TILs  are  typically  characterized  by  a  downward  phase  progression, 
which  is  similar  to  that  shown  in  Figure  15.  “Lower”  TILs  are  normally  not  correlated  in 
such  a  way,  also  in  accord  with  the  lower  range  shown  here.  It  therefore  seems  possible 
that  our  observations  can  be  explained  in  these  terms,  even  though  the  lower  range  in 
Figure  15  is  above  80  km.  (According  to  Meriwether  and  Gerrard  [2004],  the  complexity 
of  the  processes  involved  in  TIL  formation  sometimes  obscures  the  distinction  between 
the  two  types.  For  example,  variations  in  the  strength  of  the  tide  on  a  particular  day  could 
modify  the  altitude  range  in  question.)  However,  it  is  premature  to  draw  broad  conclu¬ 
sions  from  this  one  sample.  In  fact,  we  find  that  at  30°S,  TIL  altitudes  above  85  km  in  the 
early  morning — the  only  time  at  which  they  are  numerous  and  the  distribution  is 
narrow — descend  with  LT  at  a  rate  much  lower  than  that  of  the  GSWM  temperature 
maximum.  To  carry  this  further,  it  is  important  to  refine  the  process  whereby  TIL 
altitudes  are  determined.  It  is  also  important  to  develop  the  model  fully,  utilizing  com¬ 
plete  temperature  profiles  for  full  latitudinal  and  seasonal  ranges,  in  order  to  be  definitive 
about  the  information  that  SABER  yields  about  the  relation  between  the  migrating  tides 
and  TILs. 

3.  MODELING  MESOSPHERIC  BORES 

3.1.  Introduction 

Bores  are  a  class  of  waves  that  can  propagate  with  large  amplitudes  over  long  dis¬ 
tances  in  a  channel.  Their  propagation  characteristics  have  been  studied  both  analytically 
and  computationally  using,  for  example,  special  solutions  to  the  Navier-Stokes  equations 
(NSEs).  These  nonlinear  partial  differential  equations  describe  the  momentum,  contin¬ 
uity,  and  energy  relationships  that  govern  fluid  flow  in  general,  including  the  behavior  of 
air  parcels  undergoing  adiabatic  motion  in  the  case  of  atmospheric  bores.  The  influences 
of  gravity,  temperature,  and  pressure  on  the  fluid  parcels  are  also  incorporated  into  the 
NSEs.  Although  tropospheric  bores  have  been  observed,  studied,  and  modeled,  it  has 
only  been  relatively  recently  that  airglow  emissions  in  the  mesosphere  have  been  imaged 
and  associated  with  the  passage  of  ducted  bore  fields  [e.g.,  Taylor  et  al.,  1995;  Dewan 
and  Picard,  1998;  Dewan  and  Picard,  2001;  Smith  et  al.,  2003;  She  et  al.,  2004].  The 
ducting  regions  in  these  cases  are  understood  to  be  formed  by  a  strongly  stable  layer, 
above  and  below  which  are  weakly  or  neutrally  stable  regions.  Such  guiding  structures 
can  be  caused  by  temperature  inversion  layers  and  certain  types  of  wind  profiles. 

Although  there  have  been  investigations  of  bores  in  other  contexts  using  assumed 
analytical  models  for  both  the  background  wind  profile  UQ(z)  and  the  Brunt  frequency 

N(z),  we  were  interested  in  using  actual  mesospheric  measurements  of  key  physical 
variables  such  as  temperature  and  wind.  N  in  this  report  is  the  natural  buoyancy 
frequency  at  which  air  parcels  oscillate  parallel  to  a  vertically  directed  gravitational  field 
and  is  determined  by  the  temperature  and  the  lapse  rate.  The  data  used  in  this  report  were 
lidar  measurements  [ She  et  al.,  2004]  of  temperature  and  background  wind  profiles  for  a 
bore  field  being  ducted  in  the  mesosphere  and  also  included  the  observed  bore  speed 
from  a  CCD  imager.  The  data  and  the  bore  speed  observed  by  these  authors  provide  a 
basis  against  which  predictions  from  a  mathematical  model  can  be  compared.  Our  first 
objective  was  thus  to  develop  a  space-time  evolving,  ducted  mesosphere  bore  field  model 
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and  to  compare  the  predictions  from  it  with  their  observations.  Having  such  a  model,  one 
could  then  use  its  evolving  displacement  amplitude  to  quantitatively  estimate  the  strength 
of  airglow  emissions  at  any  altitude  in  the  mesosphere. 

One  approach  for  creating  the  desired  mathematical  model  is  to  pose  and  solve  an 
initial-value  boundary- value  problem  involving  the  fluid  dynamical  NSEs,  since  they 
define  the  physics  governing  the  motion  of  air  parcels.  Following  Ono  [1975],  we  used  a 
weakly  nonlinear  analysis  and  obtained  a  separation  of  variables  solution  to  the  NSEs. 
This  led  to  a  description  of  vertical  air  parcel  displacements  q(x,z,t)  caused  by  the 
passing  bore,  in  which  s(x,z,t)  is  found  to  be  proportional  to  the  product  of  the  functions 
<f>(z)  and  A(x,  t).  4>(z)  satisfies  the  Taylor-Goldstein  (TG)  equation  and /!(:/;,  t)  emerges 
from  the  solution  to  the  Benjamin-Davis-Ono-Burgers  (BDOB)  equation.  More  will  be 
said  about  both  equations  shortly. 

Ono’s  analysis  did  not  meet  our  needs  entirely;  it  lacked  the  effect  of  a  background 
wind  profile,  and  his  derivation  was  based  on  the  assumption  of  a  purely  /incompressible 
fluid.  The  effect  of  compressibility  was  needed  for  our  model  since  the  mesosphere  fluid 
is  composed  of  air  parcels.  Compressibility  was  eventually  incorporated  into  our  analysis 
using  an  approach  taken  by  Doviak  et  al.  [7997]  and  Christie  [1989]  whereby  temperature 
and  density  are  replaced  by  potential  temperature  and  potential  density,  respectively.  This 
will  also  be  discussed  shortly. 


3.2.  Formulation  of  Taylor-Goldstein  and  Helmholtz  Equations 

Pursuing  an  Ono  type  of  derivation  for  an  incompressible  fluid  [Ono,  1975],  but  now 
including  an  arbitrary  background  wind  profile,  two  basic  equations  emerge:  the  TG  and 
BDO  equations.  The  latter  is  the  same  as  the  BDOB  equation  but  without  the  Burgers 
loss  term  added  in.  Compressibility  influences  will  be  accounted  for  in  another  way. 

The  TG  equation  [Doviak  et  al.,  1991;  Christie,  1989]  eigenmodes  are  assumed  to  be 
governed  in  the  long-horizontal-wavelength  limit  by 


fpp* 

(««,') 

(t/o(z)-c0) 

-  0 ,  — oo  <  z  <  oo  . 


(1) 


(Primes  indicate  differentiation  with  respect  to  z  in  this  report.)  <j)(z)  is  the  long  wave 
eigenfunction  and  c0  is  the  eigenvalue  phase  speed  associated  with  the  bore.  Both  are  to 
be  determined.  Homogeneous  boundary  conditions  of  the  form  +  <f)  =  0  at 
z  =  -o o  and  at  z  =  oo  are  imposed  for  the  solution  of  Equation  1 .  In  practice,  we  find  that 
there  is  no  significant  difference  in  the  value  of  c0  for  the  lowest  modes,  relative  to  other 

effects  and  approximations  inherent  in  our  formulation,  whether  one  uses  Dirichlet, 
Neumann,  or  mixed  boundary  conditions. 

For  reasons  discussed  in  the  next  section,  the  TG  equation  is  transformed  to  a 
Helmholtz  equation  using  the  exponential  stratification  factor  contained  in 
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<(>  =  W  exp 


-  \dz'!H(z' ) 
V2  z 


with  H  = 


.A. 

I 

Po 


Substituting  into  the  TG  equation,  one  arrives  at  the  Helmholtz  equation 
W"  +  k2W  =  0,  with  the  propagation  constant  kz(z)  given  by 


N 2 


P0U0  M  Pq  \  +  2H' 


(U0(z)-c0f  ( U0(z)-c0 )  4  H2 


(2) 


(3) 


In  this  equation  and  following  the  lead  of  Doviak  et  al.  [1991],  the  presence  of  N2  indi¬ 
cates  that  compressibility  was  accounted  for  by  relating  density  (p0) ,  potential  tempera¬ 
ture^),  and  Brunt-Vaisala  frequency  squared  ( N 2),  respectively,  as  follows: 


— — =  -^L- ,  with  p0  — >  (const)O  1  (z) ; 

Po  6  S 

namely,  the  logarithmic  derivative  of  the  density  was  replaced  by  the  potential  temper¬ 
ature  and  its  derivative.  Figure  16  shows  the  resulting  k\  in  the  data  window  for  the  case 
of  3-point  smoothing  of  UQ(z) .  Several  positive  regions  indicate  the  presence  of  multiple 
mesosphere  ducts  that  allow  for  propagating  modes  to  be  supported. 

The  wind  profile  UQ(z )  used  in  Figure  16  and  the  calculations  discussed  herein  are 

from  the  work  of  She  et  al.  [2004].  The  data  were  originally  sampled  at  0.5  km  altitude 
increments,  subsequently  smoothed,  and  then  spline  interpolated. 


3.3.  Altitude-Limited  Data  Considerations 

Several  complicating  factors  arose  in  the  modeling  effort  associated  with  the  altitude- 
limited  data  from  the  lidar.  The  formulation  for  the  mathematical  model  requires  that  all 
terms  in  the  TG  equation  be  specified  for  the  entire  altitude  range  spanning  (0,oo)  or,  to 

an  excellent  approximation^— oo,oo).  However,  ground-based  sodium  lidar  measure¬ 
ments  of  the  event  [She  et  al.,  2004]  are  available  only  over  a  finite  altitude  range,  ~83- 
1 00  km.  Under  these  conditions,  one  can  only  hope  to  approximate  the  properties  of  the 
evolving  bore  field.  This  necessitated  making  assumptions  about  field  behavior  and 
medium  properties  outside  the  lidar  data  measurement  window. 

Another  difficulty  is  that  the  altitude-limited  measurements  of  temperature  T(z)  and 
background  wind  profile,  U0(z),  contain  noise.  Both  of  those  functions,  as  they  appear  in 

the  TG  equation,  need  to  be  differentiated  numerically  to  yield  N  ,U0  ,  and  U0  . 
Despite  initially  smoothing  both  T  and  UQ ,  data  differentiation  of  any  residual  spatial 

variations  in  these  functions  after  smoothing  was  found  to  affect  the  value  of  the  phase 
speed  eigenvalue,  significantly  so  in  some  cases  we  studied.  The  “best”  way  to  smooth 
and  differentiate  the  data  in  the  context  of  solving  the  TG  equation  has  not  been 
extensively  examined  to  the  best  of  our  knowledge.  We  experimented  with  various 
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smoothing  techniques  applied  to  the  data  of  She  et  al.  [2004],  and  settled  on  a  simple 
triangle  weighting  that  averaged  the  central  point  with  its  2,  4  or  6  nearest  neighbor 
points. 

Without  any  information  about  the  background  wind  at  altitudes  above  and  below  the 
data  window  at  the  time  of  the  event,  we  made  a  simplifying  assumption  consistent  with 
some  earlier  studies;  namely,  that  U  is  taken  to  be  spatially  uniform,  though  not  neces¬ 
sarily  of  the  same  value  above  the  data  window  as  below.  Each  of  these  two  values  is 
pegged  to  the  wind  speed  just  inside  the  nearest  border  of  the  data  window.  The  spatial 

uniformity  of  UQ  implies  that  above  and  below  the  window,  U0  and  U0  vanish.  The 
assumption  of  spatial  uniformity  of  U  is  not,  however,  absolutely  required  [Doviak  et 

al.,  1991].  The  situation  described  therein  seemed  too  exceptional  and  unlikely  to  be  of 
interest  to  us;  specifically,  it  implies  a  functional  relationship  between  potential  tempera¬ 
ture  and  the  background  wind  profile. 

Yet  another  consequence  of  having  altitude-limited  data  is  the  matter  of  how  to 
specify  TG  equation  terms  that  govern  mode  shape  outside  the  data  window.  We 
approached  this  issue  by  first  noting  that  ducted  wave  propagation  in  general  is  more 
easily  examined  using  the  Helmholtz  equation,  W"  +  k2z(z)W  =  0  when  possible.  We 

thus  transformed  the  TG  equation  to  the  Helmholtz  equation  using  the  standard  expon¬ 
ential  stratification  factor  (see  Equation  2)  and  solved  for  eigenmodes  of  the  Helmholtz 
equation.  Once  the  solutions  were  obtained,  the  resulting  field  was  transformed  back  to 
the  TG  field  variable  <f>(z) .  Following  this  approach,  we  dealt  with  the  problem  of  a  finite 
data  window  by  approximating  the  Helmholtz  mode  outside  the  data  region  with  two 
decaying  exponentials,  one  on  each  side  of  the  window.  These  exponentials  were 
connected  to  the  solution  inside  the  data  window  at  the  borders  and  made  to  vanish 
asymptotically  for  z— »  ±oo .  Outside  the  window,  terms  in  kr  (see  Equation  3)  com¬ 
posed  of  N2 ,  U 1 ,  and  U "  vanish,  leaving  the  small  negative  factor, 
- 1 /  =  -(l  +  2H')j4H2  ,  which  if  spatially  averaged  defines  a  decay  constant  in  each  of 

the  two  regions  outside  the  window.  Here,  H  is  an  altitude-dependent  pressure  scale- 
height,  the  values  of  which  were  not  available  outside  of  the  data  window  at  the  time  of 
the  bore  event  but  were  obtained,  expediently,  from  an  available  (ALOHA)  atmosphere 
model  and  some  piece-wise  temperature  fittings  using  formulas  from  the  U.S.  Standard 
Atmosphere.  While  this  would  appear  to  yield  a  gross  approximation  at  any  particular 
altitude  at  the  time  of  the  event,  as  an  averaged  quantity,  using  this  type  of  approximation 
led  to  phase  speed  results  not  too  far  from  the  observed  bore  speed,  namely,  within  1 0%. 
This  rationale  is  also  in  keeping  with  our  understanding  from  other  work  and 
observations  that  the  local  state  of  the  atmosphere  away  from  the  mesosphere  data 
window  has  a  relatively  minor  influence  on  bore  characteristics.  In  future  endeavors, 
however,  we  would  recommend  a  more  judicious  choice  of  model  atmosphere  from 
which  to  develop  H.  The  values  assigned  to  the  exponential  decay  constants  in  this 
investigation  were  thus  obtained  by  spatially  averaging  —  a2  for  each  of  the  two  regions 
outside  the  data  window. 
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3.4.  Solution  of  Taylor-Goldstein  and  Helmholtz  Equations 

Both  from  physical  considerations  and  because  the  data  used  was  space-limited,  we 
sought  a  solution  method  emphasizing  data  in  the  window  and  de-emphasizing  or 
requiring  less  information  about  conditions  outside  of  it.  To  this  end,  we  developed  a 
solution  representing  the  lowest-order  bore  modes  using  rational  Chebyshev  polynomials 
TBm(z)  [Boyd,  2000]  by  letting 


2  M 

W(z)  =  ^Cm-TBm(z).  -oo<z<co  (4) 

m=0 

Applying  the  following  cotangent-mapping,  with  #now  denoting  angle,  not  potential 
temperature, 

z--L  cot(0)  +  z  ,  with  0  <  0  <  n  ,  (5) 

and  using  a  Gauss-Lobatto  [Boyd,  2000,  p.570]  discretization  over  the  angle  via 
6  — >  0,-ln  I2M  (/ =  0,1, ...,2M),  one  finds  that  the  density  of  altitude  samples  is 
greatest  about  the  mid-altitude  z  of  the  data  window  [a,b]  =  [83,100]  km  (see  Figure  17). 
We  note  that  for  M  —  70 ,  92%  of  the  nodes  fall  within  the  data  window,  while  for  M  = 
100,  the  number  is  94%.  The  cot  # -mapping  in  this  context  also  allows  one  to  use  the 
equivalent  but  more  easily  manipulated  expansions, 

2  M  2M 

=  £c„cosm0  =>  wv,=£c„  cosmS,.  (6) 

m=0  m=  0 

The  tuning  parameter  L  can  be  chosen  in  various  ways.  For  this  investigation,  L  was 

adjusted  so  that  L= -  with  the  1st  node  being  set  at  z  =  -108A:w  «  -oo 

cot(0.5;r/2M) 

where  we  set  W  =  0. 


At  altitudes  outside  the  data  window,  the  Helmholtz  eigenmodes  were  approximated 
by  using  spatial  averages  of  u2  =  (1  +  2 H')  /  4 H2 ,  now  equal  to  a  constant  below  the 
window  and  to  another  constant  above  the  window,  as  discussed  in  the  previous  section. 
This  created  decaying  exponential  Helmholtz  mode  shapes  outside  the  data  window. 

Substituting  Equation  6  into  the  Helmholtz  equation,  after  observing  that 


d  dO  d  ,  d2 

—  - - and  — r 

dz  dz  dO  dz1 


J_  d9_ 
L  dz 


f 


1 2  > 


sin  26 - hsin2  6 — - 

dG  d62 


(7) 


yields  the  following  linear,  homogeneous  system  of  equations: 


sin2  0{ 

21} 


Lc„ 


m- 0 


w[cos(/n - 2)0,  -  cos(m  +  2)B,]  +  m 2  cos mOt 
2 

— ^-[cos(w  +  2)0  +  cos (m  -  2)0, ] 


2  M 

+k2 ( 0 \ ;c0)^Cm  cos m0,  =0  for  /  =  0, 1, .. .,  2 M. 

m=0 


(8) 
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The  method  used  to  determine  a  specific  mode  is  described  next. 

The  determinant  of  the  system  was  evaluated  numerically  for  any  specified  trial  value 
of  c0.  Stepping  through  a  sequence  of  trial  values,  and  calculating  the  corresponding 

determinant  for  each  resulting  set  of  equations,  the  c0  interval  in  which  a  root  (eigen¬ 
value)  was  located  could  be  identified  by  a  sign  change  in  consecutive  determinant 
evaluations.  A  root-finder  was  then  used  to  converge  to  the  final  phase  speed  eigenvalue. 
Having  obtained  c0 ,  the  expansion  coefficients  were  determined  in  the  standard  way  by 
normalizing  the  coefficients,  dropping  one  equation,  and  solving  the  resulting  inhomo¬ 
geneous  system.  By  repeating  this  solution  method  several  times  with  successively 
increasing  values  of  M,  a  2-pass  application  of  Aiken’s  convergence  acceleration  method 
[Press  et  al.,  1986]  on  the  rational  Chebyshev  polynomial  partial  sums  yields  phase  speed 
estimates  of  67-68  m/s.  Figure  17  shows  an  example  of  the  Helmholtz  eigenmode  thus 
obtained.  The  convergence  of  the  predicted  phase  speed  versus  M  is  shown  in  Figure  18. 
Also  shown  in  Figure  18  is  the  result  of  a  curve  fitting  experiment  using  tank  functions  to 
express  cQ(M)  at  successively  increasing  values  of  M.  Extrapolating  M— >  oo  yields  the 
phase  speed  estimate  indicated  in  Figure  18,  and  tends  to  the  same  value  that  the  Aiken’s 
extrapolation  produced. 

Based  on  observations  and  measurements  by  She  et  al.  [2004],  we  were  especially 
interested  in  bores  traveling  at  or  near  74  m/sec.  We  bracketed  our  search  for  a  mode 
traveling  between  60-120  m/sec.  This  led  to  the  2nd  lowest-order  varicose-like  mode. 
Multiple  eigenvalues  corresponding  to  higher-order  modes  with  slower  phase  speeds 
were  also  detected  when  we  applied  the  same  algorithm  to  lower  speed  ranges.  In  addi¬ 
tion,  there  was  a  heretofore  unobserved  lowest-order  mode  that  appeared  at  -168  m/sec. 

Having  obtained  the  Helmholtz  mode  of  interest  and  the  associated  phase  speed  esti¬ 
mate,  the  corresponding  TG  mode  is  obtained  using  Equation  2.  Here,  however,  we  must 
account  in  some  manner  for  losses  due  to  the  effects  of  molecular  viscosity,  eddy 
dissipation,  and  thermal  conduction  above  the  data  window.  We  used  a  small,  short  range 
loss  factor  that  varies  parabolically  with  altitude  from  and  above  the  top  of  the  data 
window  for  this  purpose.  Qualitatively,  the  resulting  0(z)  appearance  now  makes  more 
sense  physically  than  it  does  without  the  loss  function. 

In  summary,  we  developed  a  new  solution  method  for  the  Helmholtz  and  TG  equation 
that  addressed  some  of  the  issues  associated  with  altitude-limited  data  resulting  from  lidar 
mesosphere  measurements.  Using  data  from  She  et  al.  [2004],  the  phase  speed  estimate 
we  obtained  with  different  wind  profile  smoothing  was  in  the  range  of  67-68  m/s  when 
one  executes  two  passes  of  Aiken’s  convergence  acceleration  scheme  on  the  partial  sum 
expansion  for  the  eigenmode  shown  in  Equation  6.  This  procedure  yields  phase  speed 
estimates  within  10%  of  the  observed  bore  speed,  which  was  -74  m/s.  The  mode  shape 
bears  a  resemblance  to  varicose  modes  observed  in  other  bore  contexts. 

3.5.  The  Benjamin-Davis-Ono-Burgers  (BDOB)  Equation  and  its  Solution 

The  BDO  equation  governs  the  horizontally  evolving  dimension/ass  stream  function 
amplitude  quantity  A(x,  t)  necessary  to  describe  the  bore  displacement.  When  loss  due  to 
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eddy  diffusivity  is  taken  into  account,  an  added  (Burgers)  term  appears,  leaving  the 
BDOB  equation  in  dimensionless  form  as 


d  A 
dT 


dA 

+  .4  —  + 

dX 


d2 

dX2 


H(A)  -  e 


d2A 
dX 2 


=  0, 


(9) 


where  H(A)  is  the  Hilbert  transform  of  A,  e  =  n  /  (3 ,  is  the  eddy  diffusivity  coeffi¬ 
cient,  and  (3  is  defined  below. 


To  solve  this  nonlinear  partial  integro-dififerential  equation,  we  embedded  a  pseudo- 
spectral  method  for  the  spatial  variation  of  A  into  a  Crank-Nicholson  implicit  time¬ 
stepping  scheme,  [Carnahan  et  al.,  1969] .  Our  choice  for  the  particular  algorithm  used  for 
the  pseudospectral  stage  was  one  developed  by  James  and  Weideman  [1992].  It  is  based 
on  a  Fourier  series/domain  truncation  methodology.  The  truncated  domain  spanned  a 
segment  of  the  duct  axis  taken  to  be  parallel  to  the  earth’s  surface.  The  Fourier  series 
were  efficiently  evaluated  using  FFT  operations. 

The  dimension/ess  variables  (A,  X,T)  can  be  related  to  dimensions/  variables 
(A,  £  =  x  —  cQt,  t )  as  follows: 


A  =  a(b  —  a)A  /  f3, 

X  =  |/(6 -a),  and  (10) 

T  =  pt  /  (b  —  a)2 . 


(b  -  a)  is  usually  taken  as  the  duct  width.  For  several  adjacent  ducts,  as  shown  in  Figure 
16,  the  definition  of  “duct  width”  becomes  somewhat  obscure.  For  the  data  window 
altitude  range  [83,100]  km,  we  used  a  value  of  17  km.  We  note  that  |  =  x  —  cQt  is  the 

position  variable  in  the  Rayleigh  frame  moving  at  the  phase  speed  c0  with  respect  to  the 
laboratory  frame. 

In  the  above  formulas,  the  parameters  a  and  (3  arise  in  the  derivation  for  the  BDOB 
equation  and  are  given  by  the  following  formulas:  a  -a2!  a,  and  /?  =  a3  /  a, ,  where 


a , 


=  jdz 


a 
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^0  co 


+- 
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(11) 
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and 


&)=  Po(b)</>2(b)  + p0(a)<p2(a)  .  (13) 

The  BDOB  equation  thus  makes  use  of  the  previously  calculated  eigen-pair 
{4>{z)  ;cQ}  through  the  a,  (3  parameters.  Although  the  above  formulas  do  not  appear  to 

resemble  the  ones  obtained  by  Doviak  et  al.  [1991],  we  are  able  to  show  that  Equations 
1 1  and  12  agree  exactly  with  their  results  for  an  arbitrary  wind  when  the  Scorer  approx¬ 
imation  is  invoked. 

It  is  useful  to  integrate  by  parts  some  of  the  integrands  shown  in  Equations  1 1  and  12. 
This  allows  one  to  a)  minimize  the  presence  of  noisy  higher-order  derivatives  in  the 
numerical  evaluation  of  the  integrals  and  b)  to  take  into  account,  in  an  approximate 
manner,  compressibility  effects  as  done  by  Doviak  et  al.  [1991]  using  gpQ ' jp^  — >  -N1 . 

After  integrating  by  parts,  one  is  left  with  boundary  terms  and  integrals.  Outside  the  duct 

and  at  the  boundaries,  those  multiplied  by  N  ,U0  ,  or  U0  must  be  set  to  zero  in  keeping 
with  the  modeling  assumption  of  a  uniform  wind  below  and  above  the  duct  as  well  as  the 
vanishing  of  N 2  outside  the  duct.  When  this  is  done,  we  verified  that  the  vertical  and 
horizontal  parcel  displacements,  as  the  bore  passes  any  given  point  in  its  path,  are 
qualitatively  and  quantitatively  reasonable. 

Figures  19-24  illustrate  A(X,T)  as  observed  in  the  Rayleigh  frame  for  several  bores 
having  different  initial  pulse  amplitudes.  For  a  small  initial  amplitude,  e.g.,y4(W,0)  =0.1, 
as  seen  from  Figures  21  and  22,  the  bore  is  predicted  to  travel  slower  than  the  long  wave 
phase  speed,  at  least  over  the  time  T-interval  considered.  This  is  manifested  by  the  main 
peak  moving  backwards  as  T  increases.  For  larger  values  of  T,  the  bore  speeds  up  and 
eventually  travels  faster  than  the  long  wave  phase  speed.  The  situation  for  larger  ampli¬ 
tudes  shows  that  the  bore  travels  faster  than  the  phase  speed,  at  least  for  T  >  5  (see 
Figures  19  and  20).  We  note  also  that  the  shapes  of  A(X,T )  generated  in  two  examples 
using  the  BDOB  code  was  verified  by  two  published  figures  from  Christie  [1989]. 

The  results  for  a  very  long  pulse  (shown  in  Figure  23)  illustrate  both  the  front  steep¬ 
ening  due  to  the  nonlinearity  and  the  rear  flattening  caused  by  dispersion  effect  arising 
from  the  Hilbert  transform  term.  These  effects  are  also  apparent  from  the  Figure  24 
snapshot,  where  one  can  observe  the  trailing  amplitude  going  negative. 

Summarizing,  we  have  obtained  solutions  to  the  BDOB  equation  using  pseudo- 
spectral  and  Crank-Nicholson  methods,  the  results  of  which  have  been  confirmed  with 
previously  published  figures  found  in  Christie  [1989],  We  find  that  for  a  small  initial 
amplitude  pulse,  the  main  peak  of  the  bore  travels  slower,  at  least  initially,  than  the  long 
wave  phase  speed  which  corresponds  to  the  speed  of  the  Rayleigh  frame.  To  the  best  of 
our  knowledge,  this  effect  has  not  been  previously  reported.  As  the  pulse  amplitude  is 
increased,  the  bore  moves  faster.  We  have  also  examined  situations  involving  short  and 
long  pulses,  as  shown  in  the  figures.  Steepening  of  the  front  and  flattening  of  the  rear  of 
the  evolving  field  are  plainly  visible  as  one  expects  from  a  bore  analysis  involving  the 
nonlinear  and  Hilbert  transform  terms,  respectively.  A  perhaps  even  more  realistic  anal¬ 
ysis  would  replace  the  initial  pulse,  which  appears  instantaneously  at  time  T  =  0  in  this 
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investigation,  with  one  due  to  a  source  that  develops  over  time.  This  remains  to  be 
investigated  and  would  provide  an  even  more  realistic  model  of  the  evolving  field. 
Having  obtained  both  0(z)  and  A(x,t),  the  vertical  parcel  displacements  can  now  be 
evaluated. 


3.6.  Vertical  Parcel  Displacements 


Following  the  approaches  of  Ono  [1975]  and  Doviak  et  al.  [1991],  using  a  stream 
function  formulation,  one  can  show  that  the  vertical  parcel  displacements  are  given  by  the 
simple  formula 


q(x,z,t) 


0 

acQ(b  -  a) 


4>A  , 


(14) 


in  the  limit  of  no  loss  (e  =  0).  Figures  25-28  show  examples  of  the  evolving  vertical 
displacement  for  an  arbitrarily  chosen  fixed  point  x  in  the  laboratory  frame  500  km  from 
the  front  of  an  initial  pulse  of  length  AX  =  20  or  Ax  =  (b  —  a)  AX  =  340  km  . 


Two  types  of  figures  are  shown  for  the  displacements;  one  is  based  on 
a  and  (3  containing  only  Scorer  approximation  terms,  and  the  other  also  includes  non- 
Scorer  approximation  terms.  In  the  Scorer  approximation,  one  uses  in  the  derivation 
leading  to  q(x,z,t),  an  approximation  to  the  TG  and  Helmholtz  equations;  namely, 


N2 

V" 

k-<0)2 

U0~c0 

0. 


(15) 


Non-Scorer  terms  arise  when  one  includes  all  of  the  terms  in  the  TG  or  Helmholtz 
equation.  From  our  study  we  find  that  there  is  a  significant  difference  in  the  predicted 
maximum  vertical  displacements,  which  occur  at  an  altitude  of  100  km,  when  using  the 
two  parameter  sets,  as  can  be  seen  from  Figures  25-28. 


3.7.  Summary  of  Results 

To  obtain  the  displacements,  our  analysis  and  the  ensuing  codes  were  based  on  a  deep 
fluid  model  using  altitude-limited  measurements  of  temperature  and  background  wind 
profiles.  An  initial-value  boundary-value  problem  was  posed  and  solved  for  this  purpose. 
A  separation  of  variables  solution  emerged  from  a  weakly  nonlinear  formulation 
following  the  methodology  of  earlier  investigators.  Our  investigation  differs  from  earlier 
work  by  including  the  effect  of  a  background  wind  profile  based  on  lidar  data  measure¬ 
ments  in  the  mesosphere.  The  separation  of  variables  approach  required  us  to  solve  the 
TG  equation  via  a  Helmholtz  equation  solution  for  the  mode  shape  and  long  wave  phase 
speed  associated  with  the  passing  bore,  which  was  ultimately  necessary  for  parcel 
displacement  calculations.  A  rational  Chebyshev  polynomial  method  was  developed  for 
the  Helmholtz  equation  because  of  its  emphasis  on  altitude  data  samples  within  the 
mesosphere  data  window.  A  rationale  for  treating  the  region  outside  the  data  window  was 
also  developed.  The  predicted  phase  speed  was  found  to  be  within  10%  of  the  observed 
bore  speed.  We  then  chose  a  pseudospectral/Crank-Nicholson  solution  method  for  sol¬ 
ving  the  BDOB  equation,  this  being  the  other  equation  emerging  from  the  separation  of 
variables  solution.  The  computed  horizontally  varying  component  of  the  bore  field  agreed 
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with  published  figures  for  test  cases.  Based  on  our  displacement  calculations,  we  find  that 
the  Scorer  approximation  might  be  too  restrictive  to  assume,  especially  if  there  are 
multiple  ducts,  as  there  appeared  to  be  in  the  data  we  used.  Pursuing  a  direct  numerical 
simulation  in  a  future  effort  would  resolve  this  uncertainty.  We  also  think  it  would  be  of 
considerable  value  to  treat  the  source(s)  of  the  bore  more  realistically  by  using  a  temper¬ 
ature  or  pressure  time-varying  source  function  rather  than  the  impulsive  source  used  in 
this  study.  Having  developed  programs  for  calculating  evolving  vertical  parcel  displace¬ 
ments  at  any  altitude  in  the  mesosphere,  we  are  now  in  a  position  to  use  the  displace¬ 
ments  for  quantifying  airglow  emissions. 
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Mesopause  Altitude  Distribution,  June  &  July  2003,  Local  Time  1-2  hours 
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Figure  1 .  SABER  mesopause  altitude  distribution  as  a  function  of  latitude  for  2 1  June-9  July  2003 .  All  data  are 
for  local  time  1-2  h. 


Nighttime  Mesopause  Altitude  Distribution,  July  5,  2003 


-i — i — i — r 


i  i  i  i  i  ■  i 


I  i  i  r  i - 1  i  i  i  i  |  i 


E 

a> 

T3 

ZJ 


0) 

(/> 

ZJ 

(0 

CL 

o 

V) 

a) 


no  h 


100  b 


90 


80  b 


00  V 

o  o  5 


°< 

o°o« 

O  °«*  ° 


70 


—I - L - 1 - 1  I  I  I 


I  I _ I— 


-40 


60 


80 


-20  0  20  40 

Geographic  Latitude  (degrees  N) 

Figure  2.  Similar  to  Figure  1,  but  using  all  ascending  data  from  5  July  2003  only.  Local  times  were  all  between 
20  h  and  2  h.  Mesopause  altitudes  calculated  from  mean  profiles  for  5-degree  latitude  bins  are  shown  in  red. 
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Nighttime  Mesopause  Altitude  Distributions,  June-July  2003 
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Figure  3.  Distribution  of  nighttime  mesopause  altitudes,  by  latitude,  in  late  June  and  early  July  2003 
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Mesopause  Altitude  Distribution,  July  2003,  Local  Time  15-16  hours 
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Figure  4.  SABER  mesopause  altitude  distribution  as  a  function  of  latitude  for  5-17  July  2003.  AJ1  data  are 
for  local  times  15-16  h. 


2685  events,  LT  15-16  hrs  ! 
Orbits  8515-8694 
July  5-17,  2003 


o 


July  3,  2003,  Orbit  8492,  Kinetic  Temperature  Profiles 


Figure  5.  Temperature  profiles  for  three  successive  scans  on  3  July  2003.  The  effect  of  wave  structure 
on  the  definition  and  interpretation  of  the  mesopause  altitude  is  apparent. 
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Mesopause  Temperature  Distribution,  June  &  July  2003,  Local  Time  1-2  hours 
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Figure  6.  SABER  mesopause  temperature  distribution  for  the  conditions  of  Figure  1. 


Nighttime  Mesopause  Temperature  Distribution,  July  5,  2003 


Figure  7.  SABER  mesopause  temperature  distribution  for  the  conditions  of  Figure  2.  Temperatures  derived 
from  mean  profiles  for  each  5-degree  latitude  bin  are  shown  in  red. 
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Nighttime  Mean  Temperature  Profiles,  July  5,  2003 
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Figure  8a.  Mean  temperature  profiles  for  ascending  data  on  5  July  2003,  5-degree  latitude  bins  35-75°N. 
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Nighttime  Mean  Temperature  Profiles,  July  5,  2003 
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Figure  8c.  As  in  Figure  8a,  but  for  latitude  bins  5-45°S. 
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Ascending-Descending  Temperature  Comparison,  5°N-5°S,  March  27,  2002 


Kinetic  Temperature  Difference  (K) 

Figure  9.  Comparison  of  mean  temperature  within  five  degrees  of  the  equator  on  27  March  2002.  Profiles  (left) 
were  obtained  by  averaging,  separately,  the  temperature  for  all  events  on  the  ascending  and  descending  sides,  for 
which  LT~0  h  and  9  h,  respectively.  Error  bars  are  the  rms  deviations  about  the  mean.  Temperature  difference 
(right)  strongly  suggests  the  influence  of  the  diurnal  tide. 
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Temperature  Profiles,  Orbit  1630--March  27,  2002 
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Figure  10.  Four  nighttime  temperature  profiles  from  27  March  2002.  TDLs  appear  at  various  altitudes  during  each 
event.  The  largest  TIL  (event  4)  exceeds  100  K  in  amplitude. 


Temperature  Profiles,  July  13,  2002-Orbit  3218 
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Figure  1 1 .  Four  nighttime  temperature  profiles  from  July  1 3,  2002,  from  consecutive  events.  TILs  appear  at  and 
slightly  above  80  km  in  each  case. 
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Figure  12.  TIL  occurrence  probabilities,  by  latitude  and  local  time,  for  two  yaw  cycles  in  2002.  In  each  case,  the 
SABER  data  acquisition  period  began  just  at  equinox.  The  TIL  amplitude  threshold  was  20  K. 


•  Binned  Tl  layer  mean  altitude  distribution  (AT>20K) 
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Figure  13.  Mean  altitude  distribution  of  TILs,  by  latitude  and  local  time,  for  the  yaw  cycles  in  Figure  12. 
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TIL  Occurrence  Probability  (%),  AT  >  20  K 
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Figure  14.  TIL  occurrence  probability,  as  in  Figure  12  but  for  one  yaw  cycle  at  solstice.  Northern-hemis¬ 
phere  summer  is  included.  No  data  are  available  for  the  areas  appearing  in  white. 


Figure  15.  Mean  TIL  altitudes,  plotted  against  local  time,  for  equatorial  events  during  March-May  2002.  Also 
shown  are  the  location  of  local  tidal  temperature  maxima  as  predicted  by  GSWM.  For  higher-altitude  inversions, 
TIL  altitudes  decrease  with  increasing  local  time,  but  lower-altitude  TILs  appear  at  approximately  the  same 
height  regardless  of  LT. 
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Figure  16.  Helmholtz  equation  propagation  constant  (squared).  Regions  with  positive  values  correspond  to  ducts. 
Background  wind  was  smoothed  with  a  3-point  triangle  taper. 
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Figure  17.  Helmholtz  equation  eigenmode  #2  and  phase  speed,  from  the  rational  Chebyshev  polynomial  solution 
(M  =  70).  Background  wind  was  smoothed  with  a  3-point  triangle  taper. 
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Phase  Speed  Estimates  vs  M 
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Figure  1 8.  Convergence  of  phase  speed  estimates  using  rational  Chebyshev  polynomials. 
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Figure  19.  Solution  A(X,T)  of  the  BDOB  equation.  Time-staggered  snapshots  are  shown  for  T  =  0(5)100  in  the 
Rayleigh  frame.  Initial  dimensionless  pulse  amplitude  =  0.5;  eddy  diffusivity  related  loss  constant  for  Burgers 
(dimensionless)  term  =  0.1.  Note  the  forward  displacement  of  the  peak  for  increasing  T. 
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BDOB  Snapshot:  T  =  100 
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Figure  20.  Snapshot  of  the  solution  A(X,  T  =  0, 100)  of  the  BDOB  equation.  Initial  dimensionless  pulse  amplitude 
=  0.5 ;  eddy  diffusivity  related  loss  constant,  £  =0.1.  Note  the  forward  movement  of  the  main  peak.  For  very  large  T 
(not  shown),  the  peak  diminishes. 
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Time-Staggered  Snapshots  -  (Dimensionless  Variables) 
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Figure  21.  Time-staggered  snapshots  of  the  solution  A(X,T)  of  the  BDOB  equation.  Times  shown  are  T  =  0(5)100, 
and  the  solution  was  obtained  in  the  Rayleigh  frame.  Initial  dimensionless  pulse  amplitude  =  0.1;  eddy  diffusivity 
related  loss  constant  for  the  Burgers  (dimensionless)  term  =  0.1.  Note  the  backward  displacement  of  the  peak  for 
increasing  T. 
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BDOB  Snapshot:  T  =  100 


Figure  22.  Snapshot  of  the  solution  A(X,  T  =  0, 100)  of  the  BDOB  equation.  Initial  dimensionless  pulse  amplitude 
=  0.1;  eddy  diffusivity  related  loss  constant,  £  =0.1.  Note  the  backwards  movement  of  the  main  peak.  For  large  T 
(not  shown),  it  slowly  turned  around  and  started  to  move  forward  again. 
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Time-Staggered  Snapshots  -  (Dimensionless  Variables) 
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Figure  23.  Solution  A(X,T)  of  the  BDOB  equation  assuming  a  very  long  initial  pulse.  Time-staggered  snapshots  are 
shown  for  T  =  0(5)100  in  the  Rayleigh  frame;  dimensionless  pulse  amplitude  =  0.5;  eddy  diffusivity  Burgers 
(dimensionless)  term  =  0.1. 
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Vertical  Displacement  £at  xQ  =  500  km 


Figure  25.  Evolving  vertical  displacements  at  various  altitudes,  from  the  Taylor-Goldstein  and  BDOB  equations.  No 
background  wind  smoothing  was  used;  only  Scorer  terms  were  retained  in  a  and  j3 . 
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Vertical  Displacement  £  at  xQ  =  500  km 


Figure  26.  Evolving  vertical  displacements  at  various  altitudes  obtained  from  the  Taylor-Goldstein  and  BDOB 
equations.  No  background  wind  smoothing  was  used;  Scorer  and  non-Scorer  terms  were  retained  in  a  and  /3 . 
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Figure  27.  Evolving  vertical  displacements  at  various  altitudes  obtained  from  the  Taylor-Goldstein  and  BDOB 
equations.  Three-point  smoothing  of  the  background  wind  was  used;  only  Scorer  terms  were  retained  in  a  and  p  . 
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Vertical  Displacement  £  at  xQ  =  500  km 


Figure  28.  Evolving  vertical  displacements  at  various  altitudes  obtained  from  the  Taylor-Goldstein  and  BDOB 
equations.  Three-point  smoothing  of  the  background  wind  was  used;  Scorer  and  non-Scorer  terms  were  retained  in 
a  and  (3 . 
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